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The avalanche properties of models that exhibit 'self-organized criticality' (SOC) are still mostly 
awaiting theoretical explanations. A recent mapping (Europhys. Lett. 53, 569) of many sandpile 
models to interface depinning is presented first, to understand how to reach the SOC ensemble and 
the differences of this ensemble with the usual depinning scenario. In order to derive the SOC 
avalanche exponents from those of the depinning critical point, a geometric description is discussed, 
of the quenched landscape in which the 'interface' measuring the integrated activity moves. It turns 
out that there are two main alternatives concerning the scaling properties of the SOC ensemble. 
These are outlined in one dimension in the light of scaling arguments and numerical simulations of 
a sandpile model which is in the quenched Edwards-Wilkinson universality class. 



PACS numbers: PACS numbers: 05.40.+j, 05.70.Ln 



I. INTRODUCTION 

The concept of reaching "self-organized" criticality 
(SOC) in a model without any apparent tuning param- 
eter draws still attraction In many real-life systems 
power-law probability distributions are met, and it is thus 
an obvious question to ask why they would resemble or- 
dinary critical phenomena. To this end, sandpile mod- 
els have been the prevalent theoretical playground for 
the last fifteen years. Only recently their understand- 
ing has finally started to take shape. Two supporting 
approaches have been developed. The crucial notion is 
that the critical state in these models draws from both 
the boundary conditions and the drive, and also has a 
generic, field-theoretical description. This can be for- 
mulated as a variant of the directed percolation-style 
Reggeon field theory (RFT) Q , or as a mapping to in- 
terface depinning ||. The gist of this mapping 
is based on the description of the history of the sandpile 
model and its dynamics via a stochastic differential equa- 
tion, the quenched Edwards- Wilkinson (qEW) equation 
[0, H . Likewise a suitable RFT for sandpiles includes by 
necessity a density conserving term that accounts for the 
effects of diffusional transport of particles or grains. 

In this article the description of sandpiles through in- 
terface depinning is reviewed. The issue of particular 
interest to us is the physics of driven interfaces that de- 
scribe sandpile models. To this end Section II introduces 
in a short fashion the mapping, and discusses the ensem- 
ble in which SOC is reached. It is seen that this is not any 
of the normal ones familiar in the context of the qEW, 
say. The next section is devoted to a discussion of the 
standard properties of interface scaling at the depinning 
transition, and, correspondingly, the scaling laws usually 
formulated for sandpiles depending on the ensemble. We 
next concentrate in particular on a geometric description 
of depinning. This is the essential issue in obtaining the 



scaling exponents of SOC avalanches: the extension of 
the theory of 'elastic pinning paths' to the SOC case. In 
section IV numerical results are considered for the one- 
dimensional qEW sandpile in the SOC ensemble. This is 
the simplest model system in which one can try to extract 
a correspondence with the avalanche and the depinning 
pictures. This is since the ID pinning paths can be dis- 
cussed without the introduction of extra, independent 
exponents. Finally, section V finishes the paper with a 
summary about SOC properties based on recent advances 
including the numerical work presented here, and about 
remaining open problems. 



II. MAPPING SANDPILES TO INTERFACES 

Consider a sandpile model with each site a; of a hyper- 
cubic lattice having z{x,t) grains. When z{x,t) exceeds 
a critical threshold Zc{x), the site is active and topples 
so (y is a nearest neighbour of x) 

z{x, t + 1) = z(x, t) - 2d, z{y, t + 1) = z{y, t) + l. (1) 

Zc{x) is taken to be a random variable, chosen from a 
probability distribution P{zc) again after each toppling. 
If there are no active sites in the system, one grain is 
added to a randomly chosen site, z{x,t) z{x,t) + 1. 
This is the SOC ensemble, defined via the drive and the 
open boundaries, characterized by avalanches that may 
ensue after a grain has been deposited. The dynamics 
of this model can be reinterpreted through an 'interface', 
or history H{x,t) which follows the memory of all the 
activity at x. H counts topplings at site x up to time t, 
and has the dynamics 



H{x,t + l) 



H(x,t) + 1, f{x,t)>Q, 
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Hix,t), f{x,t)<0, 
This can be written as a discrete interface equation 
A H 

-^^0{f{x,t)), 



(2) 



(3) 
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with d{f) the step- function forcing it so that the interface 
does not move backwards. The 'local force' is 



f{x, t) = z{x, t) - Zc{x) 



^c(x), (4) 



in terms of (grains added to site x up to time t) 
and (grains removed from x). The fluxes and 
n°"* can be worked out in terms of the local height field 
H[x,t), and a columnar force term F{x,t) which counts 
the number of grains added to site x by the external 
drive, 



f{x, t) = V^H + F{x, t) ~ z,(x, H). 



(5) 



The step- function, 9{f), in Equation (||), the condition 
that the interface does not move backwards, introduces 
an extra noise term a{x, H) - the velocity of the interface 
is either one or zero - but for the current example this 
should not be relevant. Combining the three sources of 
effective noise, F, Zc and ct, one ends at the discretized 
interface equation 



V^H + 7]{x,H) + F{x,t)+aix,H). (6) 



Here the quenched noise r/{x,H) ~ ~Zc(x,H), and 
we have obtained the central difference discretization 
of a continuum diffusion equation with quenched noise, 
called the linear interface model (LIM) or the quenched 
Edwards- Wilkinson equation ^ . The continuum limit 
reads dth{x,t) = L'\/^h{x, t) + F + r]{x, h{x, t)). 

The SOC ensemble is illustrated in Figure (|l|). Once 
the local force is increased, by adding a grain and mak- 
ing F{x,t) — > F{x,t) -I- 1 at the chosen x an avalanche 
starts since the force overcomes the pinning force, V^i?+ 
ri{x,H) + F{x,t) > 0. The interface moves at x by one 
step, AH{x) — 1. In the subsequent dynamics of the 
avalanche the columnar force term F{x) does not change. 
For the SOC sandpiles, the correct choice of the interface 
boundary condition is H = which is to be imposed at 
two "extra sites" {x = 0, x = L + 1 for a system of size 
L in ID). The ever-increasing (Fit)) leads to an aver- 
age parabolic interface shape via the cancellation of the 
Laplacian by F. It is to be stressed that this implicitly 
contains the physics of the SOC state: the driving force 
F{x,t) is increased so slowly that the avalanches do not 
overlap and are therefore well-defined B. 



III. SCALING PROPERTIES OF ENSEMBLES 

Equation (^ exhibits a depinning transition at a 
threshold force Fc in a normal ensemble. The interface 
configuration and dynamics develop critical correlations 
in the vicinity of the critical point. For the case of point- 
correlated disorder the normal way to analyze the LIM 
is to use the functional renormalization group method. 
One-loop expressions for the exponents are found in pa- 
pers by Nattermann et al., Narayan and Fisher ||, and 




in the more recent ones by Le Doussal and collaborators 
§. 

The problem is technically and from the fundamen- 
tal viewpoint difficult, since the whole disorder correla- 
tor is renormalized. The point-random field noise term 
forms one of the universality classes in this problem |j ; 
others are the LIM with columnar noise 
quenched Kardar-Parisi-Zhang equation |11 , 
in more general terms. The mapping of SOC models to 
variants of LIM holds however some surprises. For in- 
stance, while the conjecture that the Manna model jl^ 
should be in the usual point-disorder LIM class seems to 
be true in 2D, this is clearly not so in ID if considered 
in the depinning ensemble p6| . The complications arise 
since an arbitrary choice of sandpile rules can lead to 
non-standard noise correlations that do not need to have 
a priori the same RG fixed point as the random field or 
force case. 

In the interface language the relevant exponents to de- 
scribe the depinning phase transition are: v (the correla- 
tion length exponent), z (the dynamical exponent), and 
X, the roughness exponent. Moreover assuming that the 
avalanche dynamics suffices to describe the interface dy- 
namics off the critical point, 9 = v{z — x) holds for the 
velocity exponent, with v ^ {F — FcY 0, |^. For point- 
like disorder the first-loop functional RG results cited 
above read x = (4 - d)/3, and z = 2 - (4 - c;)/9. No- 
tice the exponent relation v = 1/(2 — x) which manifests 
together with the 6'-exponent relation the fact that there 
is only one temporal and one spatial scale at the critical 
point. 

The typical quantity to be measured in the interface 
context is the interface width w (mean fluctuation) which 
in most cases equals other measures like the two-point 
correlation functions and the structure factor . From 
the sandpile viewpoint, these measure the correlations 
and fluctuations in the activity history, that is in the 
avalanche series. Initially the width grows, as a func- 
tion of time, as w t^ defining the 'growth exponent' 
[3 until either saturation is reached (with a finite order 
parameter, w), or the interface gets pinned at or below 
Fc- Scaling now implies [3z = x, as for general interface 
models. 

The LIM obeys an important invariant, with the static 



response scaling as g x(9i'^ = 0) ^ g , so that 
7/i/ = 2. 



(7) 



For forces below Fc, the (bulk) response of the interface 
triggered by a small increase in F is Xbuik = d (H) / dF ~ 
(Fc — F)~^ . By assuming that the avalanche due to a 
point seed scales as AH ^s^£^,D = d + x (since H 
scales with the roughness as £^), a hyperscaling relation 
can be derived for 7. Right at the critical point |j, 
the roughness of the interface scales as £^ and assuming 
that A (H) will scale in the same way it follows 



7 = 1 + Xi^- 



(8) 
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This also implies x + l/i/ = 2, as noted above. The 
standard scaling relations are valid for parallel dynamics: 
all sites with dH/dt > are updated in parallel. For 
extremal drive criticality (updating one unstable site at 
a time) the dynamic exponent reads zed — 2. 

The ID LIM is a bit more peculiar than one might 
expect. Numerics (which has recently been matched by 
the 2-loop RG results of LeDoussal et al.) implies that 
X ^ ■ ■ ■ 1.25 which is larger than the 1-loop and scal- 
ing argument result Xid = 1- The physical interpreta- 
tion for the fact that x > 1 has been dubbed "anomalous 
scaling" |l^ , and arises from a divergent mean height 
difference between neighboring sites, with t — > oo. 

The SOC steady-state is characterized by the probabil- 
ity to have avalanches of lifetime t and size s which follow 
power-law distributions: p(t) = t~'^* ft{t / L") and p{s) = 
s-^f{s/L"), with s - t^/^ and zin - 1) = D{t - 1). 
One can also characterize the avalanche by its linear di- 
mension, p{e) = i-^<'fe{i/L), with Te = 1 + D{t - 1). 
Here the size scales as s ^ £^ and the (spatial) area as 
i"^ (for compact avalanches) with £ the linear dimension. 
The fact that each added grain will perform of the order 
of [|l9| topplings before leaving the system leads to 
the fundamental result 



is) ^ 



(9) 



independent of dimension Ig^l- Thus, r = 2 — 2/D and 
Tt = I + {D - 2)/z. Equation (||) yields ^jv — 2, where 
7 describes the divergence of the susceptibility (bulk re- 
sponse to a bulk field) near a critical point, if one gen- 
eralizes the exponent relation of the depinning ensemble 
to the SOC case. In the particular case of an increasing 
drive and a bulk dissipation, which induces a term —kH 
to the qEW equation, this should, evidently, be valid pi[. 



IV. THE ONE-DIMENSIONAL QEW 
SANDPILE: GEOMETRIC DESCRIPTION OF 
AVALANCHES 

To go beyond the scaling exponents to a description of 
the probability distributions at the critical point of any 
particular ensemble is a more challenging task. This is 
easiest in one dimension, which we thus discuss here. The 
most well-known case in which the geometry of the ran- 
dom quenched landscape allows to use a self-affine picture 
of the progress of an interface through it is given by di- 
rected percolation depinning (DPP), the quenched KPZ 
cousin of the qEW 0, [l2[ ||. 

In this language the interface moves, e.g. in the case of 
an applied extremal drive, via a succession of punctua- 
tion events. The interface is assumed to invade the voids 
of a connected network in each of these events, in 'bursts'. 
In between these events, the pinned interface is mapped 
into a connected path on the backbone of a suitable (elas- 
tic) percolation problem For DPD the analogy is 
more or less clear, and for the qEW one speaks of elastic 
pinning paths. Such paths have the characteristics that 



the RHS of the qEW is always negative semi-definite, ie. 
/ < 0. It is not known rigorously whether such paths 
at criticality follow strictly the DPD -like scaling prop- 
erties. There are two fundamental issues: the geometric 
properties of avalanches (the scaling of voids, or whether 
the relation of size vs. area can be characterized with a 
local roughness exponent xioc) and the probability to in- 
duce an avalanche, or punctuation event if the interface 
is pushed at any par ticular spot (this issue is still open 
to discussion, see |23, |4[ p6|). 

Assuming that the DPD analogy works ||2^ , it follows 
in the depinning ensemble for the avalanche size exponent 



:,dep = l + (l/(l+X;oc))(l-lH 



(10) 



which using eg. the exponent relation v = 1/(2 — x) and 
a reasonable value for x/oc, taken to equal the global x, 
produces 



1.08. 



(11) 



For the SOC ensemble, the description of the criti- 
cal state is given in terms of the avalanche distribution 
exponents - note that due to the inhomogeneity of the 
ensemble e.g. perturbing the system off the critical point 
is more complicated [||. One has then to ask the ques- 
tion what is the prediction of pinning-path picture. Due 
to the symmetries (static response) of the qEW it could 
be assumed that the typically parabolic interface profile 
is irrelevant, except perhaps for finite size corrections in- 
duced by the boundary condition that is imposed on the 
pinning paths since H ^ at the edges. This would im- 
ply, in particular, that for the SOC avalanche exponent 
it is found that 



Ts,SOC 



(12) 



A parallel approach is to use the invariant (|9|) which leads 
to the prediction 



Ts^soc ^2-2/D = 2- 2/(1 + xioc) 



(13) 



by using the roughness exponent to estimate the cut- 
off dimension of avalanches. This is identical with the 
pinning path estimate above. 

These are then the major issues: is the depinning geo- 
metric description applicable to the SOC ensemble and if 
not why? The answer should depend only on fundamen- 
tal similarities or differences between the ensembles, and 
not on the particular model - class of qEW-like models - 
nor the dimensionality. The study of the outcome is the 
easiest in ID whereas in higher dimensions further inde- 
pendent exponents are needed, i.e. assumptions to de- 
scribe the probability distributions since the avalanches 
have in addition to an area vs. volume relation also a 
perimeter length vs. area one p5[ . 

To study the issue we next outline some numerical re- 
sults on the ID qEW/LIM, obtained from a Leschhorn- 
like cellular automaton [pTl |23] for the interface problem. 
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The system is run as a sandpile in that after the inter- 
face gets pinned a new avalanche initiation is done by 
increasing the local force Fi at a random location by one. 
System sizes upto L = 8192 have been studied, with 2 
X 10^ avalanches for the largest sizes The resulting 
avalanche size distributions, after logarithmic binning, 
imply that the effective r-exponent varies with L, and 
that the weight of the power-law -like tail increases with 
the size. Meanwhile, the effective roughness exponent 
slightly decreases. Due to the systematic finite size cor- 
rections it does not make sense to try a normal datacol- 
lapse using a fixed D and tsoc for all L. This is a little 
bit surprising given that relation (|^) is fulfilled within nu- 
merical accuracy by the data, and that higher momenta 
of the size distributions indicate just simple scaling of 
the avalanche size distribution. Certainly more numeri- 
cal analysis is called for, but two facts are worth pointing 
out. First, even for the largest system size the effective r 
is way off from the predicted 1.08 (1.024 for L = 8192). 
Second, by blindly - without any a priori justification - 
using the scaling Ansatz 



Tsoc{L) = Tsoc{L = oo) -I- At(L), 



(14) 



with a best-working function At, one observes that there 
is an apparent power-law correction to the exponent. 
This extrapolates very slowly in L to tsoc(oo) ~ 1.115, 
ie. clearly off the dcpinning ensemble value. From tsoc 
it would be in principle possible to derive the other ex- 
ponents as well, using the consequences of Eq. (^) , 



V. CONCLUSIONS 

Above, we have discussed a strategy to understand so- 
called self-organized criticality by mapping the history 
of a sandpile model to a driven interface in a random 
medium. For SOC models this idea is useful since it helps 
to understand universality classes (via noise terms gen- 
erated) though there are theoretical challenges in the un- 
derstanding of the possible classes (RG fixed points), and 
in the role of the discretization (called ct- noise above). 
Such work follows the historical connections of SOC to 
depinning, and extends it by explicitly constructing the 
right ensemble to reach SOC, and by outlining a gen- 
eral strategy to understand various models. Applied to 
other absorbing state phase transitions, in general, the 
history/interface description should be of interest. In 
some cases (e.g. the contact process) it can give rise to 
new, seemingly independent exponents psf . 

Once one has defined the right ensemble for SOC in in- 
terface depinning, the most pertinent question becomes 
if the usual avalanche exponents in SOC can be derived 
from those of the depinning transition. Here we have ad- 
dressed the question, but since the outcome is still open, 
it is worth reiterating the two possibilities. Either even- 
tually, by studying large enough systems numerically, the 



exponents of the statistically homogeneous ensemble are 
recovered once finite size effects become negligible. Or, 
it becomes apparent that the SOC ensemble is an inde- 
pendent one. The microscopic reason for this would be 
that the density of grains (average force for the interface) 
is non-uniform: in one dimension it is easy to see that 
a site X will get a larger grain flux from its neighbor on 
the bulk side and a smaller one from the boundary side 
- more trivially, there is a net flux of grains towards the 
boundaries. This inhomogeneity may persist in the ther- 
modynamic limit, in which case the avalanche relations 
will be determined by its scaling properties. 

Consider the idea depicted in Figure (||). The in- 
tegral of the average force deviation at x, AF{x) — 
Fdep - F{x)soc, where Fdep and F{x)soc are averages 
at the critical points of the ensembles, can be used to de- 
fine a finite-size correction to the normal critical point, 
as J AF{x)dx = 6F{L). It is seen that 5F will be depen- 
dent on the exact scaling function of F(x)soCy comput- 
ing which is thus the crucial issue. It will give indirectly 
the true correlation length exponent vsoc in the SOC 
ensemble via the L-dependence of the finite size correc- 
tion, which may or may not be the same as for the de- 
pinning case. The implication can be rephrased so that 
the usual exponents like v are derived from the RG in 
an ensemble with statistical translational invariance in 
the z-direction. It is thus not obvious whether the prop- 
erties of self-organized critical state actually follow from 
boundary -induced criticality or are related to the usual 
depinning one. In this respect the usual SOC models 
discussed here are inherently more complicated than the 
boundary driven cases like the so-called rice-pile model. 
This is in particular true if the symmetries of the depin- 
ning transition are broken by the SOC ensemble, as is 
the case for the quenched KPZ equation JT4[ . 

To summarize, the above problem is central in un- 
derstanding SOC-like systems. It may also be tackled, 
perhaps, from the viewpoint of absorbing state phase 
transitions and their field-theoretical description which 
provides a parallel route to the interface one. There are 
many other interesting issues, like the early-stage dynam- 
ics (growth exponent /3, suitably defined for the SOC en- 
semble), the question of the pinning path/manifolds in 
higher dimensions and so forth. For all these it is invalu- 
able that one can resort to continuum descriptions of the 
SOC sandpiles, that also seem to be inter-connected, in 
an intriguing fashion |29, 30 1. 
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H (change) 

H(x) 

H (before) 




H(0) = H(L+1) = 

FIG. 1: One-dimensional schematic example of the interface 
or history representation of a SOC sandpile model. The mean- 
field interface follows an average parabolic shape, which also 
implies that in the SOC steady-state v{x) is parabolic. No- 
tice the boundary conditions H = that ensure the loss of 
particles (equalling the increased elastic energy) in that state. 



H(x) 



H (changed) 



H (before) 




avalanche, dH(x) 



H(1) = H(L) H(1) = H(L) 

FIG. 2: An avalanche at the critical point of the normal, 
translationally invariant, depinning ensemble. The statistical 
description of such avalanches follows the scaling law (s) = 
Z^^''-, where I is (here in ID) the area of the avalanche and s 
its size. 



<F(x)>, 



AF(x) = F(x)^^p-F(x)^^^ 



soc 



FIG. 3: The average force per site, {f{x)) in the dcpinning 
and SOC ensembles. (/) is to be computed by computing the 
number of grains per site between the avanlanches. 



